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Many protocols for quantum computation require a quantum memory element to store qubits. 
We discuss the accuracy with which quantum states prepared in a Josephson junction qubit can 
be stored in a nanoelectromechanical resonator and then transfered back to the junction. We find 
that the fidelity of the memory operation depends on both the junction-resonator coupling strength 
and the location of the state on the Bloch sphere. Although we specifically focus on a large-area, 
current-biased Josesphson junction phase qubit coupled to the dilatational mode of a piezoelectric 
nanoelectromechanical disk resonator, many our results will apply to other qubit-oscillator models. 
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PACS numbers: 03.67.Lx, 85.25.Cp, 85.85.+J 



I. INTRODUCTION 

Nature provides many ways to realize a single qubit, 
but quantum information processing is useful only if 
many qubits can be coupled in such a way that two-qubit 
operations can be performed. Because qubits must be co- 
herent yet controllable, the macroscopic quantum prop- 
erties and long coherence times of superconductors make 
Josephson junctions strong candidates. 1 " Several pro- 
posed quantum computing architectures involve coupling 
Josephson junction (JJ) flux, phase, or charge qubits to- 
gether with LC resonators fi*2*MiSiSi2i2i2iiSiii supercon- 
ducting cavities mechanical resonators jiLi£ 
or other types of oscillatorsii2i22i2ii2£ Such resonator- 
based coupling schemes have the advantage of additional 
functionality resulting from the ability to tune the qubits 
relative to the resonator frequency, as well as to each 
other. Although harmonic oscillators are ineffective as 
computational qubits, because the lowest pair of levels 
cannot be frequency selected by an external driving force, 
they are quite desirable as coupling elements. 

In architectures based on JJs coupled to resonators, the 
resonators store single qubit states, transfer states from 
one JJ to another, entangle two or more JJs, and mediate 
two-qubit quantum logic. In effect, the resonators are the 
quantum computational analog of the classical memory 
and bus elements. In this paper we discuss the speed and 
accuracy with which a state can be stored in a resonator 
and later retrieved, which depends on both the state be- 
ing stored and on the coupling strength between the JJ 
and the resonator. The model we consider includes no 
disssipation or decoherence, and any loss of fidelity is a 
consequence of purely coherent quantum dynamics of the 
coupled qubit-oscillator system. This issue is essential to 
consider when designing a real quantum computer. 

The specific architecture we consider is a large-area, 
current-biased Josephson junction phase qubit coupled 
to a nanoelectromechanical resonator .iLiS The phase 
qubit is attractive not only for its robust quantum co- 
herence, but also because an effective method for state 
preparation, manipulation, and measurement has been 
developed^ Our results are obtained by direct numeri- 
cal integration of the time-dependent Schrodinger equa- 



tion with a Hamiltionian that is analogous to that of a 
two-level system in an electromagnetic cavity, and many 
of our result will apply to other qubit-oscillator models. 



II. SUPERCONDUCTING PHASE QUBIT 
COUPLED TO NEMS RESONATOR 



The low-energy dynamics of a JJ is determined by 
the difference tp between the phases of the spatially uni- 
form order parameters in the superconductors forming 
the junction. The Hamiltonian for the system we con- 
sider is H = Hj+H Ies +SH, where Hj = ~E c -^ + U(ip) 
is the Hamiltonian of the J J with current bias lb, with 
U = -Ej(cos(p+ sip) and s = h/h- E c = (2e) 2 /2C is 
the charging energy and Ej = HIo/2e is the Josephson 
coupling energy, with C the junction capacitance and Iq 
the critical current. In the large- area J J of interest here, 
Ej is much larger than E c . uj p0 = ^2E c Ej /h is the 
zero-bias plasma frequency. The lowest two eigenstates, 
|0) and |1), are used to make a qubit. -ff res = Hu>oa ! a 
is the Hamiltonian for the resonator, with a 1 and a the 
creation and annihilation operators for dilatational-mode 
phonons of frequency luq ■ The resonator is a piezoelectric 
disk sandwiched between two capacitor plates. Finally, 
the interaction term is 5H — —ig(a — a')ip, where g is 
a coupling constant with dimensions of energy that de- 
pends on the geometric and material properties of the 
resonatorii£ii& This Hamiltonian, which is analogous to 
a few-level atom in a single-mode electromagnetic cavity, 
is discussed further in Refs. 01 an d [13 ■ 

The junction Hamiltonian Hj depends on the dimen- 
sionless bias current s, which is time-dependent. We ex- 
pand the state of the coupled system in a basis of instan- 
taneous eigenstates \mn) s of Hq = Hj + H res , defined by 
H (s)\mn) s = E mn (s)\mn) s , where \mn) s = \m)j®\n) les , 
with |m)j and |ro) re s the eigenstates of the uncoupled JJ 
and resonator, respectively. The wave function is then 
expanded as 



1 ' ,J *o y '\mn) 



(1) 



The probability amplitudes in the instantaneous interac- 
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tion representation satisfy 



ihc mn — ^""^ (mn\SH — ihdt\m'n') 



x e (i/h)J^dt'[B mn (a)-E m/n ,( a )] c 



(2) 

All effects of dissipation and decoherence are assumed to 
be negligible over the time scales studied here. 

We will assume that the JJ states are well approxi- 
mated by harmonic oscillator eigenfunctions, which is an 
excellent approximation unless s is very close to unity. 
Transitions between the instantaneous eigenstates are 
caused by nonadiabatic variation of s, through the term 
(mn\-7^\m'n') s = (mn\-^\m'n') s s, which can be evalu- 
ated analytically in the harmonic-oscillator limit. 18 



III. NEMS RESONATOR AS A QUANTUM 
MEMORY ELEMENT 

An arbitrary qubit state 

|^)j = a|0)j + /3|l)j (3) 

prepared in the JJ can be stored in the ground and one- 
phonon states of the resonator's dilatational mode as fol- 
lows: Assuming the resonator is initially in the ground 
state and the JJ is detuned from the resonator, the cou- 
pled system is prepared in the initial state 

(a|0)j + /3|l)j)®|0) res = a|G0) + 0\1Q), (4) 

with |a| 2 + |/3| 2 = 1. The bias current s is then adi- 
abatically varied to tune the JJ level spacing Ae = 
huj p o(l — s 2 ) 1 / 4 to fkoo, reaching the resonant value 
s* = y/1 — (cjo/^po) 4 at time t = 0. Neglecting any 
nonadiabatic corrections, the probability amplitudes in 
the instantaneous interaction representation at this time 
are 



Cmn(O) = (a 6 m + (3 S m l)S n0 . 



(5) 



If the interaction strength g is small compared with 
the level spacing Ae, the subsequent dynamics is well 
described in the rotating-wave approximation (RWA) of 
quantum optics^ In this approximation, and on reso- 
nance, we can write Eq. (J2J) as 



9 r 
n 



f Vn + l X 1 C ,n+1, 



(0) 



where xqi = (0|<p|l)j = £* /^/2 is an effective dipole mo- 
ment. Here I* = (2E C /Ej)i(l - s* 2 )~s is the character- 
istic width in ip of the harmonic oscillator eigenfunctions 
in the JJ, when tuned to the resonator. Using Eqs. (JHJ 
and jnj we obtain, for t > 0, 

coo(i) = a 



coi(t) = (3 sin(f ) 



cio(i) 



P cos(^) 



and all c mn {t) with n > 1 equal to zero. f2 = 2gxm/K is 
the resonant vacuum Rabi frequency. 

After a time At = ir/£l, the nonvanishing probability 
amplitudes are 



coo(t) 
coi (t) 



a 
P, 



(8) 



corresponding to the interaction-representation state 
|0)j <8> (a|0) rcs + /5|l) rcs ). In this sense, the qubit state 
of Eq. 10) has been stored the in the resonator's vacuum 
and one-phonon states |0) res and |l) rcs - The JJ is now 
adiabatically detuned from the resonator. 

To retrieve the stored state, we again bring the systems 
into resonance at time t\. Using the stored amplitudes 
of Eq. © as initial conditions, and assuming that t > t\, 
the RWA equations now lead to 



coo(i) 
coi(i) 
cio(i) 



= (3 cos[§(t-ti)] 
= -0am[%(t -*i)], 



(9) 



the others vanishing. This time the systems are held in 
resonance for an interval At = 37r/f2, after which the 
original state (0J is recovered. 

The above analysis, which is based on the adiabatic ap- 
proximation and the RWA, suggests that perfect memory 
performance can be obtained with arbitrarily fast gate 
times. This is incorrect, of course, because the actual 
quantum memory performance is controlled by the cor- 
rections to these approximations, which we shall analyze 
numerically in the following section. 



IV. QUANTUM MEMORY FIDELITY 

The accuracy of a storage and retrieval operation can 
be characterized by the absolute value of the overlap be- 
tween the intended and achieved final states, or the fi- 
delity F of the memory operation. Accounting for the 
fact that the intended and actual final state have the 
same phase factors resulting from the time-evolution of 
the instantaneous eigenstates, F is given by the absolute 
value of the inner product of the intended and achieved 
interaction-representation probability amplitudes. In our 
results we actually plot the fidelity squared, 



F 



I a* coo (<f) + P*c w (t { ) 



(10) 



(7) 



which characterizes the •probability that the memory de- 
vice operates correctly. The fidelity will have a weak 
dependence on tf, which we shall discuss below. 

In Fig. we show the results of simulating the stor- 
age and retrieval of the qubit state 2 — a (|0) + |1)), which 
is on the equator of the Bloch sphere. The JJ is that 
of Ref. [23|], with parameters Ej = 43.05 meV and E c — 
53.33 neV, and the resonator has a dilatational-mode fre- 
quency Uq/2tt of 15 GHz. g is O.ObhujQ. The dimension- 
less bias current on resonance is s* = 0.545, and the 
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FIG. 1: Storage and retrieval of the state 2" 2 (|0) + |1)). The 
solid curve is the overlap squared with the inital state. After 
60 ns the qubit is successfully retrieved with a squared fidelity 
of about 94%. The dotted curve gives the occupation of the 
state 2 — 2 ( |00) + |01)) in which the qubit is stored in the res- 
onator. The dimensionless interaction strength g/hwo here is 
5%. The dashed curve is the bias current s. (inset) Expanded 
view of the weak tf-dependence of the fidelity. 



off-resonant value of 0.407 was determined by optimizing 
the fidelity. The memory operation of Fig.^was achieved 
with a squared fidelity of about 94%. 

The inset to Fig. ^ shows the final stages of the oc- 
cupation of the state 2 — 2 (|00) + 1 10) ) . Unless the state 
being considered happens to be an eigenstate of the cou- 
pled JJ-resonator system, the fidelity will depend on tf. 
However, the variation is quite small, so it is reasonable 
to assign an average fidelity, which is what we do below. 



A. Memory fidelity versus coupling strength 

The memory operation shown in Fig. Stakes two vac- 
uum Rabi periods — about 45 ns — to complete, disregard- 
ing times during which the system is detuned. The gate 
speed can be increased by increasing g, but then the fi- 
delity decreases. In the upper panel of Fig.[2]we show the 
memory fidelity squared for the qubit state 2~2 (|fj) + 11)) 
as a function of g/hujQ. As expected, the fidelity gradu- 
ally decreases with increasing g. However, there are de- 
viations from strictly monotonic dependence on 5, which 
we believe to be consequences of our use of a fixed (not 
optimized for each g) nonresonant value of s. The lower 
panel of Fig. Ogives the gate time as a function of g/HujQ, 
again disregarding nonresonant evolution. These results 
suggest that memory fidelities better than 90% can be 
achieved using phase qubits and resonators with coher- 
ence times of a few hundred ns. 
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FIG. 2: (upper panel) Memory fidelity for equator state 
2~i (|0) + |1)) as a function oig/hu>o, in percent, (lower panel) 
Time needed to store and retrieve state as a function of g/huio. 



B. State dependence of memory fidelity 

Interestingly, the memory fidelity depends not only on 
the strength of the JJ-resonator interaction, but also on 
the qubit state itself. This is because, as we explained 
above, the fidelity is determined by the corrections to the 
adiabatic and rotating- wave approximations, and these 
corrections are state dependent. Using a Bloch sphere 
representation 



cos(f)|0> 



sim 



^)e^|l) 



(11) 



for the stored qubit, we show in Fig.|3|the memory fidelity 
along two great circles, from |0) to |1) along = (left) 
and around the equator (right) starting and finishing at 

2-*(|0> + |l». 

The dependence of F 2 on 9 can be understood as fol- 
lows: When 6 = 0, the initial state of the coupled sys- 
tem is 1 00), because the resonator always starts in the 
ground state. For a weakly coupled system, 1 00) is close 
to the exact ground state for any s, because there are 
no other \mn) states degenerate with |00). In the adia- 
batic ds/dt — > limit, the large component of 1 00) in the 
exact instantaneous ground state will remain there with 
unit probability, a consequence of the adiabatic theorem, 
leading to a high memory fidelity for the qubit state |0). 
The |1) state, by contrast, derives no protection from 
the adiabatic theorem and is subject to errors caused by 
the corrections to the RWA. The weaker <j> dependence, 
which favors equator states pointing in the 2~2 (|0) + |1)) 
direction, or in a spin-i language with |0) = 1 1) and 
|1) = II): equator states pointing in the x direction, is 
not understood at present. 



V. DISCUSSION 

We have explored the accuracy with which piezoelec- 
tric nanoelectromechanical resonators can be used to 
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FIG. 3: State dependence of memory fidelity, for the same JJ- 
resonator system studied in Fig. with g/hu)o = 0.05. (left) 
Fidelity squared as a function of 8, along the arc cj> — on 
the Bloch sphere, (right) Fidelity squared as a function of <j> 
around the equator. 



store qubit states prepared in a current-biased Joseph- 
son junction. We find that the memory fidelity depends 



on both the resonator- JJ interaction strength as well as 
the position of the state on the Bloch sphere. Finding 
the optimum balance between fidelity and gate opera- 
tion time will be essential in the design large-scale su- 
perconducting quantum computers. Overall, our simu- 
lations suggest that after further optimization, generic 
states on the Bloch sphere should be able to be stored 
and retrieved in a few hundred ns with accuracies better 
than 90%, which would be a significant accomplishment 
in experimental quantum information processing. 

We expect that many of the results presented here 
will apply to other qubit-oscillator systems as well. In 
particular, the memory operations discussed here could 
be carried out with the superconducting transmission 
line resonator architecture currently being developed at 
Yaleiii^ 
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